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The production of a sterile species via active-sterile mixing in a thermal medium is studied in an 
exactly solvable model. The exact time evolution of the sterile distribution function is determined 
by the dispersion relations and damping rates Fi^ for the quasiparticle modes. These depend on 
7 = r aa /2AE, with F aa the interaction rate of the active species in absence of mixing and AE 
the oscillation frequency in the medium without damping. 7 < 1,7 > 1 describe the weak and 
strong damping limits respectively. For 7 <C 1, Ti = r aa cos 2 # m ;F2 = r aa sin 2 # m where 9 m is 
the mixing angle in the medium and the sterile distribution function does not obey a simple rate 
equation. For 7 > 1, Ti = Y aa and Y2 = r aa sin 2 2# m /4-y 2 , is the sterile production rate. In this 
regime sterile production is suppressed and the oscillation frequency vanishes at an MSW resonance, 
with a breakdown of adiabaticity. These are consequences of quantum Zeno suppression. For active 
neutrinos with standard model interactions the strong damping limit is only available near an MSW 
resonance if sin 28 <C a w with 9 the vacuum mixing angle. The full set of quantum kinetic equations 

Ofor sterile production for arbitrary 7 are obtained from the quantum master equation. Cosmological 
resonant sterile neutrino production is quantum Zeno suppressed relieving potential uncertainties 
\ associated with the QCD phase transition. 
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I. INTRODUCTION 



D 

Sterile neutrinos, namely weak interaction singlets, are compelling candidates to explain a host of cosmological and 
astrophysical phenomena. They could be a suitable warm dark matter component P, [J H, 0) IE B 
£NJ . 03], may also be relevant in the latest stages of stellar collapse [HI, primordial nucleosvnthesis|17l Il8l|. and provide 
a potential explanation for the anomalous velocity distributions of pulsars [l9l I20I l2l| . Although sterile neutrinos are 
ubiquitous in extensions of the standard model[22, 23 L 24], [2f| , the MiniBooNE collaboration[26j has recently reported 



results in contradiction with those from LSND|27I. |2E| that suggested a sterile neutrino with Am 2 ~ 1 eV 2 scale. 
Although the MiniBooNE results hint at an excess of events below 475 MeV the analysis distinctly excludes two 
neutrino appearance-only from — > v e oscillations with a mass scale Am 2 ~ 1 eV 2 , perhaps ruling out a light sterile 
neutrino. However, a recent analysis [29j suggests that while (3 + 1) schemes are strongly disfavoured, (3 + 2) neutrino 
schemes provide a good fit to both the LSND and MiniBooNE data, including the low energy events, because of 
the possibility of CP violation in these schemes, although significant tension remains. These issues notwithstanding 
the MiniBooNE result does not constrain a heavier variety of sterile neutrinos such as those that could be suitable 
warm dark matter candidates with masses in the keV range [11, [U, [H, 0, H HI, H EH E3- Their radiative decay would 
contribute to the X -ray background [H 13, [3(| HH HH from which constraints on their masses and mixing angles may 
be extracted d, l3ll. l32[ [H[. It has also been suggested that precision laboratory experiments may be sensitive to 
~ keV neutrinos |34j. Being weak interaction singlets, sterile neutrinos can only be produced via their mixing with 



an active species, hence any assessment of the possibility of sterile neutrinos as dark matter candidates or their role 
in supernovae must begin with understanding their production mechanism. To be a suitable dark matter candidate, 
two important constraints must be satisfied: the correct abundance and a velocity dispersion that restricts the free 
streaming length to be consistent with the constraints from structure formation. Both ingredients depend directly on 
the distribution function of the sterile neutrinos, which in turn depend on the dynamics of production and evolution 
until freeze-out. 

Pioneering work on the non-eq uilibrium dynamics of neutrinos in a medium was cast in terms of kinetic equations 
for a flavor "matrix of densities" [35| or in terms of 2 x 2 Bloch-type equations for flavor quantum mechanical stat es l36l . 
f37l | . A general field theoretical approach to neutrino mixing and kinetics was presented in [38|, H^] (see also [25|). 
however sterile neutrino production in the early Universe is mostly studied in terms of simple phenomenological 
rate equations [J H, 0, |4l|, EH, and numerical studies 0, [H[ rely on an approximate semi-phenomenological 



'Electronic address: boyan@pitt.edu 



2 



approach [4ll |42| . A field theoretical study of the hadronic contribution to the sterile production rate near an MSW 
resonance has been reported in ref . [44j . 

Understanding the dynamics of oscillations, decoherence and damping is of fundamental and phenomen olog ical 
importance not only in neutrino cosmology but also in the dynamics of neutral meson mixing and CP violation [45l l46l. 
\at\ and axion-photon mixing in the presence of a magnetic field [25], a phenomenon whose interest has been rekindled 
by the recent results from the PVLAS collaboration [48| (see the discussion in ref. [4^] ) . As argued in[5(3] the spinorial 
nature of neutrinos is inessential to describe the dynamics of mixing and decoherence in a medium. 

Recently we reported on a study [5 1| of mixing and decoherence in a theory of mesons that provides an accurate 
description of similar phenomena for mixed neutrinos. This effective theory incorporates interactions that model the 
medium effects associated with charge and neutral currents for neutrinos and yields a picture of the dynamics which 
is remarkably general. The fermion nature of the distributions and Pauli blocking effects can be simply accounted for 
in the final result [51]. This study implemented quantum field theory methods to obtain the non-equilibrium effective 
action for the "neutrino" degrees of freedom. More recently this approach was extended to study the production of 
sterile neutrinos both from the effective action as well as from the correct quantum kinetic equations obtained directly 
from the quantum master equation (52j. The results obtained in ref.[H2| clarify a host of important aspects, such as 
the approach to equilibrium and a detailed analysis of quantum Zeno suppression when the decoherence time scale is 
shorter than the oscillation time scale, thereby confirming previous results obtained for neutrinos with standard model 
interactions in refs. [5^, ■ The study in refs. [5ll |52| relied on integrating out the bath degrees of freedom, assumed 
to remain in equilibrium, up to second order in a perturbative expansion akin to an expansion in Gj? in the standard 
model. This perturbative treatment restricted the analysis to the weak damping regime in which the decoherence 
time scale is larger than the oscillation time scale. In refs. [12, [53| it was pointed out that a strong damping regime 
featuring the opposite relation between these time scales could emerge near an MSW resonance for small vacuum 
mixing angle consistent with constraints from the X-ray background [9ll3ll l32l. |33| . 

Motivation and goals: A sound assessment of sterile neutrinos as warm dark matter candidates requires a reliable 
description of the kinetics of production and evolution towards freeze-out. Strong departure from equilibrium in the 
distribution function at freeze-out could lead to significant changes in the abundance or skewed velocity distributions 
that could affect the free streaming lengths and structure formation 55]. In this article we complement and extend a 
previous stud y 15211 on the non-equilibrium production of a sterile species via active-sterile mixing. While the previous 
study (5ll [52l . |53| focused on the weak damping limit consistently with a perturbative expansion in standard model 
interactions, this article studies an exactly solvable model that allows to explore systematically the strong damping 
case and to draw general conclusions on the production dynamics of a sterile species. 

The model incorporates all the relevant ingredients: active-sterile mixing via a mass matrix which is off-diagonal 
in the flavor basis, and the coupling of the active species to a continuum of degrees of freedom which are taken as a 
thermal bath in equilibrium and includes an index of refraction contribution which modifies the mixing angles and 
dispersion relations in the same manner as for neutrinos propagating in a medium. 

Summary of results: The exact solution of the Heisenberg equations of motion allows a complete investigation 
of the non-equilibrium dynamics of production of the sterile species in the weak and strong damping regimes and to 
analyze in detail quantum Zeno suppression. We obtain the quantum master equation and from it the complete set 
of kinetic equations that describe the production and evolution of the active and sterile distribution functions and 
coherences and reproduce the exact results. Our main results are : 

• The exact solution of the Heisenberg (-Langevin) equations of motion for one active and one sterile species yield 
two different modes of propagation in the medium corresponding to quasiparticles whose dispersion relations and 
damping rates (widths) depend on the dimcnsionless ratio 7 = T aa /2AE with r aa the active species interaction 
rate in absence of mixing, and AE the oscillation frequency in absence of damping but including the index of 
refraction in the medium. The weak and strong damping cases correspond to 7 < 1 and 7> 1 respectively. 
The exact distribution functions for the active and sterile species are obtained, their time evolution is completely 
determined by the widths of these quasiparticles and the oscillation frequency including corrections from the 
index of refraction and damping. 

• The results in the weak damping regime 7 <C 1 coincide with those obtained previously in refs.(5ll. l52l. [53j: the 
dispersion relations are akin to those of neutrinos in a medium with an index of refraction and the damping 
rates are V± = T aa cos 2 8 m ; T2 = T aa sin 2 6 m where 9 m is the mixing angle in the medium. The generalized 
active-sterile transition probability obtained from expectation values of Heisenberg operators in the full quan- 
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turn density matrix is sin ™ m e~ Flt + e~ T2t - 2e"^( r i+ r 2)* cos [AEt] 
cannot be described by a simple rate equation, since the distribution function depends on the time scales 

i/ri,i/r 3) i/AE. 



• In the strong damping regime 7 » 1 the oscillation frequency vanishes at an MSW resonance signaling a 
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breakdown of adiabaticity, and the widths of the quasiparticles become Ti ~ T aa , T2 ~ r aa sin 2 29 m /4rf 2 . 
To leading order in 1/7, the time evolution of the sterile distribution function simplifies into a rate equation, 
with the production rate given by T 2 ~ sin 2 29 m (AE) 2 /T aa (see eqn. (|III.17|) 'I. The active-sterile transition 
probability is strongly suppressed ~ I/7 2 . The vanishing of the oscillation frequency, the suppression of the 
transition probability and the production of the sterile species are all manifestations of the quantum Zeno effect 
emerging in the strong damping limit. 

• For active neutrinos with standard model interactions it is shown that the strong damping limit is only available 
near an MSW resonance for small vacuum mixing angle 9 satisfying the condition sin 29 < a w where a w . This 
condition is likely satisfied by the constraints on the vacuum mixing angle from the X-ray background 0, [30j 
and entails that sterile neutrino production is strongly suppressed by the quantum Zeno effect near an MSW 
resonance. This suppression may relieve uncertainties from the QCD phase transition for keV sterile neutrinos. 

• The quantum master equation for the reduced density matrix is obtained under standard approximations. From 
it the generalized transition probability and the complete set of kinetic equations are obtained valid in all regimes 
of damping. These reproduce the results obtained from the exact treatment. Under simple approximations the 
full set of kinetic equations is presented in the form of quantum kinetic equations for a "polarization vector" . 
The complete set of kinetic equations (jIV.22IIIV.24p along with the relations pV.25HV.26[) provide a complete 
description of the non-equilibrium evolution of the active and sterile distributions and coherences. 



II. THE MODEL 



The main ingredients in the dynamics of the production of a sterile species via active sterile mixing are: i) a mass 
matrix off diagonal in the flavor basis which mixes the sterile and active species, ii) the coupling of the active species 
to a bath in equilibrium. In the standard model the bath degrees of freedom are quarks, leptons or hadrons, these 
equilibrate via strong or electromagnetic interactions, hence can be taken to be in thermal equilibrium. 

We propose a simple exactly solvable model that includes all these ingredients, it is a generalization of a model 
for quantum Brownian motio n l5rl |57| which has long served as a paradigm for the study of quantum dissipative 
systems in condensed matter |58| and quantum optics 59]. It consists of a set of coordinates q that describe the 
"system" coupled to a continuum of harmonic oscillators Q p that describe a thermal bath in equilibrium. This simple 
model is generalized so that the coordinates q ajS stand for the active and sterile "neutrinos" , these are mixed by off 
diagonal elements in a frequency matrix but only the "active" coordinate couples to the bath degrees of freedom. The 
motivation for studying this model stem from the realization that the spinorial degrees of freedom are not relevant 
to describe the non-equilibrium dynamics [50], a statement confirmed by previous studies of mixing, oscillations and 
decoherence in a theory mesons [U HH which yields a remarkably robust picture of the dynamics of neutrinos. 

The Lagrangian for this model is 



L 



q-q 1 fc 2 I + M 2 +V ) q 



It. 



Ql - w$Ql 



la CpQp 



where the flavor vector is given by 



q = 




and k is a momentum label, which is assumed but not included as an argument of q ay 
2x2 identity matrix and 



M 2 M 2 

aa as 

M 2 Ml 



V aa (k) 




(III) 



(II.2) 



for compact notation, I is the 



(II.3) 



The off diagonal elements of the mass matrix M lead to active-sterile mixing and the matrix V models a "momentum 
dependent matter potential" for the active species. 

A sum over k makes explicit the field theoretical nature of the model, however just as in the case of neutrinos, we 
are interested on the dynamics of a given k mode in interaction with the "bath" degrees of freedom. 

The correspondence with neutrinos is manifest by assuming that the matter potential is obtained from one-loop 
charged and neutral current contributions of 0{Gf) from a background of leptons, quarks or hadrons (or neutrinos 
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in equilibrium) and features a CP-odd term proportional to the lepton and baryon asymmetries and a CP-even 
term that only depends on energy and temperature [60l [6l| . The linear coupling of the active species to the bath 
degrees of freedom with C p oc Gf models the charged current interaction, for example the coupling between the 
electron neutrino and protons, neutrons and electrons in a medium, GfiPp(Cv — CaJ ^^tpN'ipe'JiiO- — l b )v e (see 
a similar description in(25l. l39j). The label p will be taken to describe a continuum when the density of states is 
introduced below. Obviously the model pi.l[) affords an exact solution and yields a remarkably general description 
of the dynamics. The main ingredient is the coupling of a degree of freedom to a continuum of bath or environmental 
degrees of freedom. Such coupling to a continuum is also at the heart of particle- antiparticle oscillations in neutral 

meson systems (K° — K :B° — B ) as described in refs.Bg, 47J. Other versions of this model, without mixing have 
been studied with focus on the dynamics of equilibration[62l. l63j| . 

For vanishing matter potential V the flavor q a s and the mass coordinates gx,2 are related by an orthogonal trans- 
formation 



U{6) 




U(9) 



cos 9 sin t 
— sin 9 cos ( 



(II.4) 



where the orthogonal matrix U{9) diagonalizes the mass matrix M 2 , namely 



U- 1 (9)M 2 U(9)= ( M * ° 



Mf 



(II.5) 



and 9 is the "vacuum" mixing angle in absence of the "matter potential" V. 

In the flavor basis the mass matrix M can be written in terms of the vacuum mixing angle 9 and the eigenvalues of 
the mass matrix as 



„ _ 2 SM 2 
1 = M 1 + 



-cos 29 sin26> 
sin 29 cos 29 



(11.6) 



where we introduced 



M 2 = - (Mf + Ml) ; SM 2 = Mf - M\ 



(II.7) 



The frequencies of the flavor modes are determined by the diagonal entries of the matrix M 2 in the flavor basis, 
introducing 



ZJ(k) 



k 2 + M 



(II. 



these are given by 



Lo a (k) = uj(k) 



5M 2 
2U(ky 



cos 29 



co s (k) = ui(k) 



SM 2 
2U(k)' 



cos26> 



(II.9) 



Focusing on the relevant case of ultrarelativistic neutrinos, we anticipate that the only approximation to be invoked 
is the one in which U(k) is larger than any other energy scale. It is convenient to introduce 



k 2 l- 



■V= l^fc) 2 + -^II- 



SM 2 



cos2 9--^ 
sin 29 



sin 29 
cos20- _ 



(11.10) 



The exact solution will be presented in the Heisenberg picture, in which the density matrix is time independent 
and determined by its initial value, which is assumed to be uncorrelated and of the form 



/3(0) = p q ®p Q 



(11.11) 
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The bath is taken to be in thermal equilibrium with density matrix pq = Tre~ HQ / T where Hq is the Hamiltonian for 
the sum of free harmonic oscillators of frequencies W p . 

The Heisenberg equations of motion for the coordinates q a , s , Q p are the following 

q a + K aP qp = i] a ; a,/3 = a,s (11.12) 
Qp + W^Q P = q a C p , (11.13) 

where we have introduced the flavor vector 



V : 



= J2°pQp (l) ■ ( IL14 ) 



The solution of eqn pi,13|) is 



where 



a " 



w pJ0 



p (t) = Q p °\t) + g- sin [W p (t - t')]q a (t')dt' , (11.15) 



Qf(t) = ^=[i p e-^ + 4e^] , (11.16) 

is a solution of the homogeneous equation and A p , A p are free field annihilation and creation operators with the usual 
canonical commutation relations. The distribution function for the bath degrees of freedom is 

Trp Q AlA p = — 1 = n(W p ) (11.17) 

e t — 1 

Introducing the solution (|II.15|) into (|II.12[) we find the Heisenberg-Langevin equations [59j 

q a (t)+K aP qp(t)+ f E a p{t - t')qp{t') = £ Q (i) (11.18) 



where the self energy is diagonal in the flavor basis and given by 



E«0(i -*') = - E ^ sin ^ - *')] (Jo)- ( IL1 9) 



The stochastic quantum noise is 

at) = Y,c P Q P 0) (t) (l), (n.20) 



and we note that 

Trp((t)=0. (11.21) 

The self energy E is written in dispersive form by passing to a continuum description of the bath degrees of freedom, 
writing 



C 2 

where the density of states 



C 2 f°° d 

£ -f sm[W p (t - t')} = i / — ImE M ( W )e-( t - t ') (11.22) 



7rC 2 

Im£ QQ (w) = £ — £ - W p ) - <S(^ + Wp)] (11.23) 



2W V 
p ^ 



has the properties 

ImE M (-w) = -ImE oa (w) ; ImS aa (o;) > for > . (11.24) 
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The density of states Iml] aa contains all of the relevant information of the bath. The Heisenberg-Langevin equation 
pi,18[) is solved by Laplace transform, introduce 



(11.25) 



q a (s) = / e st q a {t)dt ; etc, 
Jo 

in terms of which the equation of motion (|II.18|) becomes an algebraic equation 



s 2 S af 3 + K af: -j + £ a/ s(s) 



qp (s) = q a (0) + sq a (0) + ( a (s) 



where in the flavor basis 



In what follows we need the analytic continuation of the self-energy to real frequencies s — > ilu + + 

E QQ (s = ioj + + ) = ReS oa (o;) + HmE aa (uj) 

with the dispersive relation 



ReE aa (u;) = ~V J 



and P stands for the principal part. 

The solution of eqn. (|II.18|) in real time is given by 



with 



q a (t) = G afi (t) qf3 (0) + G a0 {t)q p {O) + / G a p{t')ip{t - t')dt' 

Jo 



ds 



G a p{t) = / — G aP (s) e st . 
Jc 27r * 



The Laplace transform of the propagator is given by 



G(s) 



s 2 I + K + S(s) 



(11.26) 



(11.27) 



(11.28) 



(11.29) 



and C is the Bromwich contour that runs parallel to the imaginary axis and to the right of all the sing 
in the complex s-plane. It follows from eqns. (|II.32|) and (|II.26[) that the propagator matrix G a p(t) is a 
solution of the equation of motion (|II.18|) with initial conditions 

G aP (0) = ; G Q/3 (0) = 1. 
It is convenient to introduce the following combinations 

1 



A( S ) 
P(s) 



SM 2 



Saa(s) + V a , 



and the matrix 



1 



(cos 26»- A(s)) 2 + sin 2 26 



cos26>-A(s) -sin 28 
-sin 29 -cos26» + A(s) 



(11.30) 



(11.31) 



(11.32) 



ularities of G 
homogeneous 



(11.33) 

(11.34) 
(11.35) 

(11.36) 
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in terms of which we find 



G(s) = 



1 



I + A(s) 



I - A(s) 



2 S 2 +W 2 {k) + ifi(A(s)-p(s)) 2 s 2 +ZJ 2 (fc) + iMi( A(s) + ~ (s) ) 



(11.37) 



Each term in this expression features poles in the complex s-plane near s ~ ± i uj{k) which are found by first performing 
the analytic continuation s — *■ iu> + + upon which the denominators in G{s) become 



s 2 + uJ 2 (k) + 5 ^f. (A(s) T P(s)) -> -uj 2 + uJ 2 (k) + ~ 



ReE aa (w) + iIm£ aa (o;) + V aa 



6M 2 
T — p(^) 



where the analytic continuations are given by 



p(w) 



(cos26»- A fl (o;) - iA/(w)) 2 + (sm26>) 2 



ReE aa (w) + y aa 



SM 2 



Aj(a;) 



ImS] aa (o)) 
<5M 2 



(11.38) 



(11.39) 



(11.40) 



The complex poles describe quasiparticles, the real part determines their dispersion relation and the imaginary part 
their damping rate in the medium. At this stage it is convenient to introduce the following variables 



A R = A R {U{k)) 



V aa + Re£ aa (cJ(fc)) 



5M 2 



(11.41) 



and write 



where 



Po 



= Aj (57(A)) = ImS aa (57(A)) 
po SM 2 p 

P {u5(k))=p Q re- la 



( cos 26» - Ar\ 2 + ( sin 26>) ' 

(l- 7 2 ) 2 + (27COs2(? m y 
27 cos 29 



a = — arctg 
2 6 



1 -7 2 



(11.42) 
(11.43) 

(11.44) 

(11.45) 
(11.46) 



and the branch is chosen such that < arctg [• • •] < 7r. The mixing angle in the medium, 9 m , is defined by the 
relations 



cos26>-A fi „ sin 20 
cos 29 m — ; sm 29 m — 



Po 



po 



an MSW resonance in the medium occurs whenever [22J, |23|, [24J 

cos 29 = Ar . 

The only approximations to be used are the following 



SM 2 RcE aa H _ ImS aa (a;) 



57(A) 



a;(A) 



o7(A) 



(11.47) 



(11.48) 



(11.49) 



these are all consistent with the ultrarelativistic limit, small radiative corrections and the narrow width limit, all 
approximations used in the case of neutrinos. Using these approximations we find the following complex poles: 
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The first term in (|II.37|) features complex poles at 



lj = ± Qi + i— - 



with 



ReEaa(^(fc)) + V aa ~ 5M 2 po T COS Ct 



Ti = 



r sin a 



(11.50) 

(11.51) 
(11.52) 



The second term in (III.37P features complex poles at 



with 



Q 2 = U(k) + 



i 



4u;(fc) 
r sin a 



to = ± + i— 



Re£ aa (w(fc)) + Vaa + 3M 2 po r cos a 



7 



(11.53) 

(11.54) 
(11.55) 



where 



ImS aa (^(fc)) 
aJ(fc) 



(11.56) 



is the interaction rate for the active species in absence of mixing in the limit ZJ(k) SM 2 , which is of relevance 
for ultrarelativistic or nearly degenerate neutrinos. In what follows we suppress the argument uJ(k) in the quantities 
Ar.i, etc., to simplify notation. 

Near the complex poles the analytic continuation G(s = ioj + + ) features a Breit-Wigner form, and the inverse 
Laplace transform can be performed by approximating the analytic continuation by the Breit-Wigner Lorentzian. We 
find 



G(t) 



e 2*1 



2if2i 2 

JQ 2 t e -^t i 

2 



2iil 2 



I + T 



2i?l[ 2 
e - l a 2 t g-^t x 

2 



I + T 



I- T* 



(11.57) 



where we have neglected wave function renormalization (residues at the poles) and introduced the complex matrix 



T = 



cos 28 m — vy 
— sin 20m 



— sin 29 rn 
cos 29 m + ij 



(11.58) 



where all quantities are evaluated at u> = u7(fc) and used the approximations (|II.49[) . Inserting the result (|II.57|) 
into the solution (|II.30j) we obtain the complete solution for the time evolution of the Heisenberg operators. The 
Breit-Wigner approximation leading to exponential damping in (|II.57[) is a Markovian approximation[59]. The full 
solution requires the initial conditions on the Heisenberg operators g(0),q(0), it is convenient to expand these in a 
basis of creation and annihilation operators of flavor states, 



ff /J (0) = -=[o^(0) + oJ(0)l ; q {O) 



/2uTr. 



a (O) - 4(0) 



(3 — a,s 



(11.59) 
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where u> aiS are the frequencies associated with flavor eigenstates given by eqn. (|II.9|) . Under the validity of the 
approximations (|II.49|1 . we can approximate 



leading to a simplified form 



q a (t) 



I + T* 



T* 



0/3(0) + /i.e. 



a/3 



Under the same approximations, we find the Heisenberg annihilation operators at an arbitrary time from 

Pa(t)~ 



these are given by 



a a (t) 



<?o(*) - 



I + T* 



; Pa(t) = q a (t), 



e -m 2 t e -^t 1 



I-T* 



O/j(0) 



a/3 



ZU(fc) 



u(fc) 



&(*-0 



a/3 



(11.60) 



(11.61) 



(11.62) 



(11.63) 



where we have used the initial condition G a p(0) = (see eqn. pi.33[l ) and (|II.49[I . Under these same approximations 
we find 



G(t 



iG(t') i 



uj(k) w(k) 



I + T* 



I-T* 



(11.64) 



A. Transition probability 

The result (jll.63j) allows us to obtain the generalized transition probability from expectation values of these operators 
in the initial density matrix. Denoting (a(t)) = Tr pa(t) and using the result (|II.21jl we find 



(a a (t)} = I e 



-tOii g-^t I 



I + T* 



e -iu 2 t e --£t_ 



,£2,1 

2 



I-T* 



(0/9(0)). 



(11.65) 



a/3 



Consider an initial density matrix that yields an initial non- vanishing expectation value for the annihilation operator 
of the active component, but a vanishing expectation value for the sterile one, namely 

(a«(0)) ^0 ; (a s (0)) =0 (11.66) 

From the form of the matrix T given by eqn. (|II.58|) we find the generalized active-sterile transition probability 



7>»-.(*) 



M*)) 



MO)) 



sin 2 26 T , 
4r 2 



e -rit + e -r 2 t _ 2e -i(r 1+ r 2 )t cos „ ^ 



where 



SM 2 p r Im£ aa (ZZ7(fc)) 
— "i = — cos a ; li+l2 = - 



2tU(fc) 



(11.67) 



(11.68) 



and Ti 2 are given by eqns. pi.52III.55H . The expression pi.67|) is similar to the transition probability for particle- 
antiparticle mixing of neutral mesons 45, 46]. 

The oscillatory term is a result of the coherent interference between the quasiparticle states in the medium and its 
exponential suppression in pi.67[) identifies the decoherence time scale r^ec = 2/(Ii + T%) = 2/T aa . 



10 



B. Weak and strong damping: quantum Zeno suppression 

The above expressions for the propagation frequencies and damping rates of the quasiparticle excitations in the 
medium lead to two different cases: 

7 <C 1 =>• weak damping (11.69) 
7 > 1 =>■ strong damping (11.70) 

These conditions can be written in a more illuminating manner, from the definitions (|II.42|) and (|II.56jl it follows that 

'-sS (,L71) 

where 

is the oscillation frequency in the medium in absence of damping, namely AE is given by |f22 — ^i| setting A/ = 0, 
i.e, the difference in the propagation frequencies only arising from the index of refraction in the medium. The 
dimcnsionless quantity 7 is the ratio between the oscillation time scale 1/ AE and the decoherence time scale 2/T aa . 
When 7> 1 the environment induced decoherence occurs on time scales much shorter than the oscillation scale and 
active-sterile oscillations are strongly suppressed. In the opposite limit 7 <§; 1 there are many oscillations before the 
environment induces decoherence. 

The strong damping condition (|II.70|) is then recognized with the condition for quantum Zeno suppression by 
scattering in a medium[4j, [36j. It corresponds to the limit in which the active mean free path is shorter than the 
oscillation length and decoherence by the medium suppresses active-sterile oscillations. 



1. Weak damping case: 7 <C 1 

For weak damping it follows that 

r«l ; sina«7COs20 TO (11.73) 
and the widths I\ 2 given by pi.52III.55p become 

rx = r aa cos 2 e m ■ r 2 = r aa sin 2 e m . (11.74) 

For the oscillation frequency we obtain 

^ = AE = W (IL75) 



and 

T 



( XX ~-ZZZ ) = ( J -1 ) »<M P.™) 



where U(8) is the unitary matrix given by eqn. (III.4j) . Introducing the Heisenberg annihilation and creation operators 
in the medium as 



(11.77) 



°i(*) ^ _ rr-ifn \ ( a a (t) \ 
a 2 (t) )- U {6m) {a s (t) ) 

and similarly with the creation operators, the time evolution (|II.65|1 in the weakly damped case yields 

M*)) ^ _ ( e- 4fil * e"^ 4 \ / (o o (0)) \ m . 

(a 2 (t)) )-{ o e -n 2f e -% t I I ( Oa (0)) ) ■ [LL ^> 
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Therefore, in the weak damping regime, the Heisenberg operators a\ 2 , Q4.2 create and annihilate the in-medium states 
that propagate with frequencies fii^ and their ensemble averages damp out with the widths I\ 2 - The active-stcrilc 
transition probability in this limit, is obtained from eqn. (|II.67j) . and is given by 



V a ^s(t) 



(a s )(t) 



K)(0) 



.-r 2 t 



2e -i(r 1+ r 2)tcos lAm] 



(11.79) 



In the weakly damped case the decoherence time scale Tdec = 2/T aa is much larger than the oscillation time scale 
\/ AE, hence many oscillations take place before the interaction with the environment leads to decoherence. 

These results reproduce those of references (5li l52l [53| and confirm their generality and applicability to the case of 
neutrinos with standard model interactions studied in ref. [53j . 



2. Strong damping case: 7 2> 1 



The case of (very) strong damping yields the following simplifications: 

r 2 ~ 7 2 - 1 + 2 cos 2 26„ 



r sin a 



sin 2 26L 



leading to the damping rates 



r 2 



27 2 



26 V, 



47 2 



sin 



2 29 n 



47 2 



(11.80) 
(11.81) 

(11.82) 
(11.83) 



This is a remarkable result, the quasiparticle width r 2 becomes vanishingly small in the strong damping regime, with 
important consequences for production of the sterile species as seen below. Furthermore, the oscillation frequency is 
found to be 



Q 2 — O x = 



5M 2 po 
2ZJ(k) 



cos2#„, = AEcos29 ri 



(11.84) 



this is another remarkable result in the strong damping regime: the oscillation frequency vanishes at the MSW 
resonance. It follows from eqns. (|II.46|) and (|II.51III.54|) that the vanishing of the oscillation frequency at an MSW 
resonance is an exact result for any 7 2 > 1. This result implies that there is a degeneracy right at the resonance, and 
unlike the quantum mechanical case in which there is no level crossing, in presence of strong environmental damping 
the two propagating states in the medium become degenerate at the resonance leading to a breakdown of adiabaticity. 
Furthermore, in this regime the transition probability pi.67jl is strongly suppressed by the factor I/7 2 -C 1, it is given 
by 



M*)) 



Mo)) 



sin 2 2B„ 

47 2 



-T 1 t 



+ e- F2t - 2e-^ ri+r2 )* cos [(Q 2 - n x )t] 



(11.85) 



In the strong damping limit Afl = |02 — Oi ] < AE hence it follows that Tdec "C 1/Att and the interference term is 
strongly damped out before one oscillation takes place. This is the quantum Zeno effect in which the rapid scattering 
in the medium prevents the build up of coherence [36j. 

The vanishing of the oscillation frequency, the suppression of the transition probability and T 2 in the strong damping 
case are all manifestations of the quantum Zeno effect. Of particular importance is the vanishing of the oscillation 
frequency at the MSW resonance because this entails a breakdown of adiabaticity. 



III. PRODUCTION OF THE STERILE SPECIES 



The number of sterile particles is given 

N s (t) = (ot(t)o.(i)) 



(HIT) 



12 



where the Heisenberg operators are given by eqns. (|II.63III.64|) and the expectation value is in the density matrix 
pi.lip . Let us consider the case in which the initial density matrix p q is diagonal in the flavor basis with initial 
populations 



N a (0) = <o+ (0)o a (0)) ; Ar s (0) = (at(0)a s (Q)>. 



(III.2) 



Using the results (|II.63|II.64|) and the stochastic noise given by eqn. (III.20III.16P with the averages (|II.17|II.21l ) we 
find 



N s (t) = V a ^ s (t)N a (0) + V s ^s(t)N s (0) + N*(t) 
where V a ^s(t) is the active-sterile transition probability given by eqn. (|II.67j) . and 



,-ifiit „-^rt 



f- + e 



-tn 2 t -±#t 



/-+ 



1 / e~ ia 
f ± = -(l±——(cos29 m + vf) 



(III.3) 

(III.4) 
(III.5) 



The contribution N^(t) is completely determined by the correlation function of the noise in the initial density matrix, 
it is given by 

sin 2 26» m f dw lmT, aa (uj) 2 
' n{uj) Fi{u;t) - F 2 (u);t) 



N s(.t) = 



4r 2 



7T 2w{k) 



where n(u>) = \e^l T — 1] 1 and 



Fi(u;t) 



LUi — OJ - 



i = 1,2. 



(III.6) 



(III.7) 



The frequency integral is carried out by approximating the functions Fi{io;t) as Breit-Wigner Lorentzians near their 
complex poles, the result is found to be 



NHt) = 



sin 2 9m cos 2 i 



^n(^)(l- e - rit ) + ^^2)(l 

M 1 2 



-r 2 t 



-e 2 



i ( r 1+ r 2 )t r aa [n(fli) + n(n 2 )] 



(l - cos(n 2 ~ Qi)tj + (fi 2 - ^i) sm{Q 2 - Qi)tj I (III.6 



The set of equations piI.3III.671 IIII.4[) and pil.8[) completely determine the time evolution of the sterile distribution 
function N s (t). 



A. Weak and strong damping limits 

1. Weak damping: I7I <C 1 



In the weak damping limit the results above yield 



ri 



1 



sin a 



0(7) ; cos a ~ 1 



, cos 9 m ; T 2 — r aa sin 



2 - Hi 



SM 2 po 



AE 



(111.9) 



2U(k) 

which lead to the following expression for the number density of the sterile species, valid for an initial density matrix 
diagonal in the flavor basis and with N s (0) = , N a (0) 7^ 0, 



N s (t) = 7V a (0) 



sin 2 26 „ 



-rit 



f e 



-r 2 t _ 26-^*008 [AEt] 

' 2, -> -»(n 2 )(i 



■ sin 2 9 m n(fii) (1 - e~ L lt ) + cos 2 9 m n(0 2 ) (1 - e~ i2t ) + 0(j) (111.10) 



this result reproduces those in ref. [52| for N s (0) — 0. We note that the production of the sterile species cannot be 
described in terms of a simple rate equation in the weak damping case because it depends on several different time 
scales. 
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2. Strong damping limit: 7 3> 1 



While the strong damping limit I7 > 1 must be studied numerically, progress can be made in the very strong 
damping regime I7I 3> 1. It will be seen below that this regime is relevant for sterile neutrinos near an MSW 
resonance. In this regime the above results yield 



r 2 r 


- 7 2 






r 2 r 


sin : 

- J-aa 4 




8M 2 p 


2ti(k) 



cos 26 m = AE cos 20 m . (111.11) 



The coefficients 



1 ri, 22 

22 = ^ — <#; 1 

2r2 /r M \ 2 x / n^ 2 7 2 + cos220 m 78 



1 r^^-fli) 2cos2(9 m 1 



< 1 (111.12) 



therefore the second line in the noise contribution ((III. 8)1 becomes subleading. Furthermore the ratios 

~ 2. < 1 

rii 7 Z 

r 4 

4^- — t — ( IIL13 ) 

r 2 T 2 sin 2 26» m 

therefore, only the term with r aa /r 2 survives in the first line in (|III.8|) . Since the transition probability in the first 
term in eqn. (|III.3|) V a ^ s cx I/7 2 (see eqn. (|II.67[1 ) this term is also strongly suppressed, therefore in the strong 
damping limit 

N s (t) ~ iVf(t) ~ n(fi a )(l - e~ T2t ) • (HI.14) 

Hence in this limit the sterile population obeys a simple rate equation 

dJV a (t ) 
dt 

however the sterile production rate is 



= T 2 [n{n 2 )-N s (t)], (111.15) 



r 2 = r aa 2^i«r aa an.ie) 
47^ 

becoming vanishingly small in the strong damping case. We conclude that sterile species production is strongly 
suppressed in the strong damping case as a conseguence of the quantum Zeno effect . The non-perturbative nature of 
this result is manifest by writing 

r 2 = sin 2 20 m -!^}!. (IH.17) 

1 aa 

We note that with 7 = T aa /2AE (see eqn. (III.71|0 this result coincides with the effective rate in the quantum Zeno 
limit 2AE/T aa <C 1 obtained in reference [42l| and implemented in the numerical study in refs. 0, [43| . However, we 
argue below that in the case of sterile neutrinos, the strong damping limit is only available near an MSW resonance, 
and far away from this resonance the non-equilibrium dynamics corresponds to weak damping and the time evolution 
of N s (t) cannot be described by a simple rate equation. 
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IV. QUANTUM MASTER AND KINETIC EQUATIONS 



Although we have obtained the time evolution of the distribution function from the exact solution of the Heisenberg- 
Langevin equations (under the approximation (|II.60[0 . within the cosmological setting it is more convenient to obtain 
a set of quantum kinetic equations for the distribution functions. This is achieved by obtaining first the quantum 
master equation for the time evolution of the reduced density matrix. In the case of neutrinos, the index of refraction 
term V aa is of first order in Gp (Fermi's effective weak coupling) while the self-energy E = S^j + iEj is of second order. 
Furthermore the study in the previous sections clearly shows that the contribution of the real part of the self-energy 
yields a second order renormalization of the index of refraction which can be simply absorbed into a redefinition of 
V aa - The most important aspect of the second order self-energy correction arise from its imaginary part, which yields 
the damping rates of the collective quasiparticle excitations. The production of the sterile species is associated with 
this imaginary part, and not the real part of the self-energy, which only renormalizes the index of refraction in the 
medium. Therefore it is convenient to include the index of refraction in the "non-interacting" part of the Hamiltonian 
by first diagonalizing the Hamiltonian for the system's degrees of freedom q corresponding to the first term in the 
Lagrangian (|II.1[> . This is achieved by introducing the mass eigenstates in the medium with the index of refraction 
as follows. The matrix K in eqn. (|II.10[) can be written as 



M 



SM 2 p 




- cos 29 m sin 29 „ 
sin 29 m cos 29 r , 



(IV.l) 



where the expressions for p and the mixing angle in the medium are the same as (|II.44[ III.47p ) but neglecting the 
second order correction ReE aa to the index of refraction. The diagonalization of the Hamiltonian is achieved via the 
unitary transformation pi.4[) but in terms of the mixing angle in the medium 9 m that includes the correction from 
the index of refraction, namely 




; U{9) 




(IV.2) 



Again to avoid proliferation of indices we refer to the coordinates that diagonalize the Hamiltonian with the index 
of refraction with the labels 1,2, which now should not be identified with those labeling the complex poles in section 

Expanding qi j2 and their canonical momenta p\^ 2 in terms of Heisenberg annihilation and creation operators 



1 



where the frequencies in the medium are 



lU 2 



di + a 



57(fc) 



Uj, 



(IV.3) 



V„„ 



uJ(k) + 



SM 2 po 
4ZU(fc) ~ 4aJ(fc) 
6M 2 p a 



(IV.4) 



4u7(fc) iU(k) 

Under the approximation (|II.49[) the active and sterile annihilation (and creation) operators a QjS are related to a\,i as 

a a = cos9 m a\ + sin0 m a2 ; a s — cos9 m a2 — sm9 m ai . (IV. 5) 

The total system-bath Hamiltonian becomes H = Hq + Hj where 



H = a l a * ^ + [ P p + W pQI 
Hi = (gi cos 9 m + q 2 sin 9 m ) C P Q P . 



(IV.6) 
(IV. 7) 



The density matrix in the interaction picture of Hq is 



p t (t) = e lHot e- lHt p{0)e lHt e~ lHot 



(IV.: 
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where p(0) is given by eqn. (|II.11[) . The equation of motion of the density matrix in the interaction picture is 

dpi (t) 



dt 



= -i[fl/(t),ft(f)] (IV.9) 



with Hj(t) — e lHot Hje lHat is the interaction Hamiltonian in the interaction picture of H n . Iteration of this equation 
up to second order in the interaction yields [5 9j| 

^ - -i [#,(*), MO)] - J dt 1 [Hj(t), [H^t'lUt')]] + ■■■ (IV. 10) 

The reduced density matrix for the system's variables q is obtained from the total density matrix by tracing over the 
bath degrees of freedom Q p , which are assumed to remain in equilibrium[59]. The following standard approximations 
are invoked[59]: a): factorization: the total density matrix is assumed to factorize 

ft(*)=p,,i(t)®Po(0) (IV.ll) 

where it is assumed that the bath remains in equilibrium, b): Markovian approximation: the memory of the 
evolution is neglected and in the double commutator in (|IV.10[) f>i(t') is replaced by f>i(t) and taken out of the 
integral [59(. Taking the trace over the bath degrees of freedom yields the quantum master equation for the reduced 
density matrix, 



dpujt) 
dt 



= - f dt'Tr PQ { [#,(*), [JT / (t'),fc(t)]] } + • • • (IV.12) 
Jo 



where the first term has vanished because TtqPq(0)Qp (t) = since Q P °\t) is a free harmonic oscillator in the 
interaction picture of Hq (see eqn. pi. 161) ). The trace over Q in the double commutator requires the following 
ingredients 

E / C A ww^ PQm ^ m ^ {t,) = E^[( 1 +"(^)) e ^ (t " t ' ) +^)^ p(t "° 

p,p' V p p' p p 

— Im£ »(l + n(w)) e -M*-f ) (IV.13) 

7T 



E 



CpCpi 



p.p 



1 V^w p w p , 



Tvp Q (0)Q p o ,\t')Q p o \t) =J^2W ( 1 + n ( W p))e~ lWpit '~ t) + n ( W p) elWAt ' 



-t) 



ImI] aa (a ) )n( W ) e -^( t -*') 



(IV. 14) 



where the interaction picture operators Q(°\t) are given by eqn. 16[) and we have used eqns. (|II.23III.24I) . 

Several standard approximations are invoked: terms that feature rapidly varying phases of the form ajaj e^i+^i) 1 

and aiOLje -1 ^ 1 ^^ 1 are averaged out in time leading to their cancellation, in the quantum optics literature this is known 
as the "rotating wave approximation" 59] , similar terms are discarded in the kinetic approach in ref. (25l . [3f| . The time 
integrals are evaluated in the Weisskopf-Wigner approximation[52, 59]. Finally we also invoke the ultrarelativistic 
approximation uj\ ~ u>2 ~ u7(fc). Neglecting the second order energy shift (see eqn. (111.29jl ). the final result for the 
quantum master equation is given by 



dpnjt) 
dt 



Y aa i 
-7TS cos 2 6 m £ii [p R ] + sin 2 9 m C 2 2 [pr] + ^ sin 



£i2[pr] + £21 [/Of 



(IV. 15) 



where Cij[pR\ are the Lindblad operators [5£ 



A?L°.r] = (l + n(u)i)) PRa\a 3 + a\aip R ~ dipRa] - ajp R a\ 
+n(u>i) pRdia] + a 3 a\p R - a\p R a 3 - a]p R a,i 



(IV.16) 
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In these expressions, the annihilation and creation operators carry the time dependence in the interaction picture, 
namely 



(IV.17) 



The trace of the reduced density matrix is automatically conserved in time as a consequence of unitary time evolution 
of the full density matrix. Denoting the expectation value of any interaction picture operator A(t) in the reduced 
density matrix by 



(A}(t)=Trp R (t)A(t), 

we obtain the following equations for the expectation values of the annihilation operators 



di 



(ai)(t) \ 

<«2>(t) J 



-iu>\ — ^f 2 - COS 2 6 r , 
sin20 TO 



-Efa sin26> m 
r • 2 
-iu>2 f 2 - sin i 



( ai )(t) 
<oa)(t) 



(IV.18) 



(IV.19) 



The eigenvalues of the matrix in eqn. pV.19[) arc found to be — zf2i,2 — Ti^/2 where fii^ are obtained from eqns. 
pi.51III.54|) by setting the second order contribution to the energy shift ReS aQ = 0, and I\ 2 are precisely given by 
eqns. pi.52III.55P but again setting Re£ aa = in p , which of course is a consequence of having neglected the second 
order energy shifts (real part of the self energy) in the quantum master equation. It is a straightforward exercise to 
obtain the (complex) eigenvectors of the matrix pV.19[) and to write (a a , s ) in terms of these through the relation 
pV.5p . Fixing the initial values of the corresponding eigenvectors to yield the initial values (a Q )(0) ^ 0; (a s )(0) = 
we find 



(a s )(t) 



(O«>(0) 



sin 2 29 r , 
4r 2 



-rit 



-r 2 t 



2e~ 



i(ri+r 2 )t 



cos 



(0 2 -Qi)t 



(IV. 20) 



which is the same as the transition probability (|II.67|) but neglecting the second order correction from Re£ aa . These 
results clearly show that the quantum master equation PV.15P correctly describes the non-equilibrium dynamics 
including the strong damping regime, the only difference with the exact result being that the second order e nerg y shift 
Rc£ aa is neglected . The quantum master equation PV.15[) is exactly the same as the one obtained in ref . [52j . 
We now introduce the distribution functions 



Uij = Tr p R (t)aj(t)aj(t) , 



(IV.21) 



the diagonal components describe the population of the in medium states, and the off-diagonal components the 
coherences [59]. Accounting for the free field time dependence of the operators a* , a in the interaction picture, we find 
the following kinetic equations for the distribution functions 



-r„ 



"ii 



h 2 2 = -T a 



"-12 



sin 



(n 2 2 - n(uj 2 )) + 



i W2-W1 mi2 



- aa 



n\2 



sin 26 „ 
4~ 

sin 28„ 
4~ 

sin 29 m 



{n\2 + n 



\n\2 



(nn + 



)} 
)} 

n 2 2 - n(cvi) - a 



12) 



l l2) 



(^2)) 



(IV.22) 
(IV.23) 

(IV.24) 



where n(u)i) are the equilibrium distribution functions. In terms of the nij(t) we obtain the time evolution of the 
active and sterile distribution functions via the relation pV.5[) , namely 



N a (t) = cos 2 9 m nii(t) + sm 2 6 m n22(t) + -sm28 m (n 12 (t)+n* 12 (t)) 



N s (t) = sm 2 9 m n 11 (t)+coa 2 e m n 2 2(t)~ -sm29 m (n 1 2{t) + n* 12 (t)) 



(IV.25) 
(IV.26) 



The weak damping limit can be studied in a perturbative expansion in 7 -C 1 by considering the terms ni2,n* 2 in 
equations pV.22|IV.23j) and the terms riu — n{uj t )\i — 1,2 in equation pV.24j) as perturbations. This study was 
carried out in ref.[52( and reproduces the result eqn. (|III.10|) for the sterile population. Therefore the set of quantum 
kinetic equations pV.22HIV.24[) reproduce the exact results both in the weak and strong damping cases. 



17 



We can now establish a correspondence with the quantum kinetic equation often quoted in the literature [25|, |42j, 
l43l . l64l |65| by introducing the following "polarization vector" [66] 



P (t) = (djjCla + a l a s)(i) = N a (t) + N s (t) 

P x {t) = {a) a a s + a\a a ){t) 

P y (t) = -i(ala s - a+a a >(t) 

P z (i) = (oio -aJa,)(t)=JV„(t)-A?.(t) 



(IV.27) 
(IV.28) 
(IV.29) 
(IV.30) 



where the creation and annihilation operators for the active and sterile fields are related to those that create and 
annihilate the propagating modes in the medium 1,2 by eqn. pV.5[) . and the angular brackets denote expectation 
values in the reduced density matrix pn which obeys the quantum master equation (|IV. 15|) . In terms of the population 
and coherences the elements of the polarization vector are given by 



Pa = nn + n 2 2 

P x = -sin20 m fnn - ri22j + cos20 m (ni 2 + n*, x 
P y = -i(nn - n* 12 ) 

P z = cos 26 m (nn - + sin 26 m (ni2 + n* l2 ) 
Using the quantum kinetic equations pV.22IIIV.24p we find 



dPo 
dt 



r r 

A aa p x aa 

~2~ z ~ ~2~ 



{nil - n(u)ifj + (n 22 - n{u 2 )^j + cos 29 m (n{wi) - n(ui 2 fj 



(IV.31) 
(IV.32) 
(IV.33) 
(IV.34) 



(IV.35) 



dP x 
dt 



-i(ui 2 - w 1 )cos20 m (n 12 - n[ 2 )) - -^P x - sin 20„ l (n(u> 1 ) - n(uj 2 jj 



(IV.36) 



dP V 
—± = -(W2 - wi)(n 12 + % 2 j - ~Y P v 



(IV.37) 



dP z ., \ • nn I * \ ^aa n Taa 

—j- = -i[u>2 - wij sin 20 m (rai2 - n 12 j —P z — 



{nn - n{ui)j + (n 22 - n(oj 2 ) 



Under the approximation oji ~ uj 2 ~ oj(fc) we can take 

(n(uji) - n(u 2 )j ~ , 

and neglect the last terms in eqns. (|IV.35|iV.36jl . Introducing the vector V with components 

V = (oj 2 — uji) ( sin 29 m , 0, — cos 29„ 
we find the following equations of motion for the polarization vector 



dP t7 - T aa f a dP „ 

— = U x P - — [P x x + P„y ) + — z . 



dt 



dt 



This equation is exactly of the form 



dP - - - dP 
— = V x P - DP T + —^z 
dt dt 



often used in the literature [36j, |42j, |43j, |6J, [65(, where 

r 

J- a. a. 



D = -^- ; P T =[P x x + P y y 



(IV.38) 



(IV.39) 



(IV.40) 



(IV.41) 



(IV.42) 



(IV.43) 
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Therefore the quantum kinetic equation for the polarization vector pV.4ip is equivalent to the full set of quantum 
kinetic equations pV.22IHV.24l) . 

However it must be highlighted that the set of equations (jiV.4lHV.42l) is not closed because it must input the time 
evolution of Po which is obtained from the full set of kinetic equations (jIV.22IIIV.24p . 

Often the last term in pV.42j) (Po) is omitted, however, such omission is not warranted, since it follows from the 
definition of P , eqn. pV.3ip and eqns pV.25|IV.26p . that 

Po = N a (t) + N s (t) , (IV.44) 

therefore Po vanishes only when both the active and the sterile species have reached equilibrium. Thus we advocate 
that the set of kinetic equations pV.22HiV.24p combined with the relations pV.25IIV.26p provide a complete description 
of active and sterile production. 



V. CONSEQUENCES FOR COSMOLOGICAL PRODUCTION OF STERILE NEUTRINOS. 

The results obtained above can be straightforwardly adapted to the case of neutrinos by replacing the equilibrium 
distributions ft(fii,2) by the Fermi-Dirac distributions in the ultrarelativistic limit and the matter potential from 
forward scattering in the medium. 

While in general 7, r li2 and fii 2 depend on the details of the interactions, masses and vacuum mixing angles, 
an assessment of the consequences of the results obtained above on cosmological sterile neutrino production can be 
obtained for an active neutrino with standard model interactions. In this case the matter potential for temperatures 
features a CP-odd contribution proportional to the lepton and baryon asymmetries, and a CP-even contribution that 
depends solely on momentum and temperature. In the ultrarelativistic limit with uJ(k) ~ k the matter potential for 
neutrinos is given by [6(| [6ll. [67|. 



L-A 



M 2 
m w 



(V.l) 



where L is proportional to the lepton and baryon asymmetries and A ~ I0[60l [6l| . for antineutrinos L — > —L. 
The active neutrino interaction rate (neglecting contributions from the lepton and baryon asymmetries) is given by 

[Ills 511, IS [MI 

T aa ~ G 2 F T 4 k . (V.2) 

For keV sterile neutrinos an MSW resonance is available only for L Th/Myy when the first term in the bracket 
in (jV.ip dominates [1, 0, |43|, [6l[, while no resonance is available when the second term dominates. We will analyze 
separately the two different cases 

L « W- (V.3) 

lvl w 

L>> & ( v - 4 ) 

1V1 W 

where we have taken k ~ T. In the first case no MSW resonance is possible for keV sterile neutrinos, whereas such 
resonance is possible in the second case [3, H, 0, |6l|. 

• High temperature limit: At high temperature above the MSW resonance for V aa 3> SM 2 and neglecting the 
second order correction to the matter potential (ReS) , 



»~Wf < v - 5 » 



For L < T 2 /M^ 

SM 2 po G F T 5 



l{k) M, 



(V.6) 

w 



and the ratio 
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where a w is the standard model "fine structure constant". For L ^> T 2 /M^ a similar analysis yields 



7 ~ G F M 2 



( T 2 



\ LM? 



w 



T 2 



< 1 



• Low temperature limit: In the low temperature regime for V aa *C SM 2 , po ~ 1 and 7 becomes 



(V.8) 



(V.9) 



however in perturbation theory V aa ^> ImE since V aa is of 0(Gp) and ImS aa ~ 0{G 2 F ). Therefore since in this 
regime 



Im£ Qa 




ImS„„ 


SM 2 p 




SM 2 



SM 2 » V aa > ImE a a 



ImE„ 



SM 2 



< 1 



(V.10) 



The conclusion of this analysis is that far away from an MSW resonance, either in the high or low temperature 
limit damping is weak, namely at high or low temperature away from the MSW resonance 



7: 



2AE 



< 1. 



Therefore the strong damping condition may only be fulfilled near an MSW resonance 
po ~ I sin20|. 



(V.ll) 
7r/4 in which case 



• Near an MSW resonance: As mentioned above a resonance is only possible for keV sterile neutrinos for 
L T 2 /M^JH, 0, [43], |61| . For very small vacuum mixing angle sin 28 <§C 1 it proves illuminating to write 
the resonance condition cos 29 — V aa /SM 2 as V aa ~ SM 2 and p ~ | sin20|, with V aa given by eqn. (|V.1|) for 
L > T 2 /M^. Therefore 5M 2 /k - G F T 3 L, hence using eqn. (|V.2|) near the MSW resonance, the ratio 



(V.12) 



8M 2 „ 



GpM 2 , 




Otyj 


( T 2 \ 


|sin26»| 




~ | sin 26» | 





Therefore, the strong damping condition near the resonance is fulfilled provided that |sin26>| -c ct w . With 
a w ~ 10 -2 the region near an MSW resonance is generally described by the strong damping regime for | sin20| < 
10 -3 , which is likely to be the case for sterile neutrinos [1, [2l| and is consistent with constraints from the X-ray 
background @, ©, [3Q, M, HI • 

In the resonance region the sterile production rate is described by the simple rate equation (see eqn. (|III.15[) ) 

N s (t) = -T 2 [N s (t) - n eq \ (V.13) 
where the sterile production rate T 2 is given by eqn. (|III.16p which can be written as 

T 2 - sin 2 20- 



(5M 2 ) 2 



w(k) 2 T aa 



(V.14) 



and clearly exhibits the suppression for small vacuum mixing angle and the non-perturbative nature as a function 

ofr QQ . 



This analysis leads to the conclusion that away from an MSW resonance the weak damping condition holds, 
sterile neutrino production cannot be described by a simple rate equation but involves Ti : 2 and AE. In this regime 
the quantum kinetic equations pV.22HV.24p may be simplified [5 2j by neglecting the terms with n\2,n\ 2 m eqns. 
PV.22IIV.23P and the terms with nu — n(u>i); 7122 — n(u>2) in eqn. (|IV.24[) . The resulting equations are very simple 
and their solutions feature the two damping rates Ti = r aQ cos 2 m ; T 2 = Taasin 2 ^. This simplification also holds 
if the lepton asymmetry is of the same order of the baryon asy mmetry L ~ 1CV 9 in which case L <§; T 2 /M 2 for 
T > 3 MeV[lI[67[ and no MSW resonance is available [43l [60L [6l| . Near an MSW resonance for sterile neutrinos with 
~ keV mass and sin 28 < 1CP 3 the strong damping condition holds and N s (t) obeys a simple rate equation, but the 
sterile production rate is suppressed by the quantum Zeno effect. 

For keV sterile neutrinos with small mixing angle sin26> < 1CV 3 , the MSW resonance occurs near the scale of 
the QCD phase transition T ~ 180 MeV[4j, |44| with the inherent uncertainties arising from strong interactions 
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and the rapid change in the effective number of relativistic degrees of freedom in a regime in which hadronization 
becomes important. However, as argued above, near the MSW resonance the strong damping condition is fulfilled 
and quantum Zeno suppression hinders the production of sterile neutrinos. As discussed above the sterile distribution 
function obeys a simple rate equation with a production rate given by eqn. (|III.16[) or alternatively 1 T[) which is 
strongly suppressed by the factor I/7 2 ~ sin 2 29 /a^ <C 1. This suppression of the sterile production rate makes the 
production mechanism less efficient near the resonance, thus relieving the uncertainties associated with the strong 
interactions, although these remain in the non-resonant scenario 



VI. CONCLUSIONS 



The production of a sterile species via active-sterile mixing has been studied in a simple, exactly solvable model 
that includes all the relevant ingredients: active-sterile mixing via an off-diagonal mass matrix and the coupling of 
the active species to a bath in thermal equilibrium. The exact solution of the Heisenberg -Langevin equations allows 
to obtain the exact time evolution of the distribution function for the sterile species and the active-sterile transition 
probability. Both are determined by the dispersion relations and damping rates (widths) of the two quasiparticle 
modes in the medium. These depend on 

7=-^ (VI.l) 
' 2AE v ' 

where r aa is the interaction rate of the active species in the absence of mixing and AE is the oscillation frequency with 
corrections from forward scattering (the index of refraction) but no damping. 7 -C 1;7» 1 correspond to the weak 
and strong damping regimes respectively. In the weak damping case the damping rates are I\ = r aa cos 2 # m ;F2 = 
r aa sin 2 9 m the active-sterile transition probability is given by eqn. (|II.79[) , and the time evolution of the sterile 
distribution function is given by eqn. (|III.10p for vanishing initial sterile population, both feature these two scales 
along with the oscillation time scale. As a result, the time evolution of the sterile distribution function does not obey 
a simple rate equation. These results confirm those of refs.[5ll [52|, The exact solution allows the systematic 
exploration of the strong damping case for which 7 > 1 corresponding to the situation in which the interaction rate 
in the medium is faster than the oscillation time scale and the quantum Zeno effect is present 36] . In this regime we 
find that the damping rates of the quasiparticles are T\ = T aa ; T2 — T aa sin 2 29 m /4rj 2 where 9 m is the mixing angle 
in the medium. The active-sterile (generalized) transition probability is 



'a — >s 



sin 2 29 r< 



4j2 



e-rxt + e -r 3 t _ 2e -^+ r ^ cos[(J}i - fi a )i 



In the strong damping regime the oscillation frequency f2i — 0,2 oc cos 29 m vanishes at an MSW resonance and the 
two quasiparticle states become degenerate leading to a breakdown of adiabaticity. The sterile distribution function 
obeys a simple rate equation with a sterile production rate T2 strongly suppressed for 7 2 3> 1. The suppression of 
the active-sterile transition probability and the sterile production rate, and the vanishing of the oscillation frequency 
in the strong damping limit are all consequences of quantum Zeno suppression. The quantum master equation for 
the reduced density matrix is derived and shown to be valid in both limits. From it we obtain the complete set of 
quantum kinetic equations that yield the non-equilibrium evolution of the active and sterile distribution functions. 
The complete non-equilibrium time evolution of the active and sterile distribution functions and the coherences are 
given by the set of equations (|IV.22HIV.24[) along with the identifications PV.25IIV.26[) . The set of kinetic equations 
PV.22IIIV.24I) are shown to be equivalent to the kinetic equations for the "polarization vector" often quoted in the 
literature. However, unlike these the set (jIV.22tiIV.24|) along with pV.25IIV.26[) yield a complete description of the 
non-equilibrium dynamics amenable to a straightforward numerical analysis, the extrapolation to fermionic degrees of 
freedom is a straightforward replacement of the equilibrium distribution functions by the Fermi-Dirac distributions. 
Furthermore, the analysis based on the exact solution and the quantum master equation yield a wealth of information 
that cannot be easily gleaned from the set of kinetic equations, for example the active-sterile transition probability. 

For active neutrinos with standard model interactions it is shown that the weak damping limit describes the 
parameter range away from an MSW resonance and that the strong damping limit only emerges near the resonance 
for very small vacuum mixing angle, such that sin 29 < a w ~ 10~ 2 . Such small value is consistent with constraints 
from the X-ray background. This result bears important consequences for cosmological sterile neutrino production. 
In the resonant production mechanism of ref.[l| the production rate peaks at the MSW resonance, however our 
analysis, which includes consistently the damping corrections, shows that quantum Zeno suppression hinders the 
sterile production rate near the resonance. For keV sterile neutrinos the MSW resonance occurs in a temperature 
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range too close to the QCD phase transition. Hadronization and strong interactions lead to substantial uncertainties 
during this temperature regime which translate into uncertainties in the production rate. Quantum Zeno suppression 
of the production rate in this regime relieves these uncertainties. 

In summary: The set of kinetic equations pV.22HIV.24p (with Fermi-Dirac equilibrium distributions) along with 
the relations (|IV.25IIV.26[) yield a complete description of the non-equilibrium dynamics of active and sterile neutrino 
production valid in the weak and strong damping limits. Quantum Zeno suppression is operative near an MSW 
resonance and suppresses the sterile production rate, thus relieving potential uncertainties associated with the QCD 
phase transition for keV neutrinos. 
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